arXiv: 1503.02223vl [cond-mat.dis-nn] 7 Mar 2015 


Statistical features of magnetic noise in mixed-type impact fracture 

Zs. Danku^, Gy. B. Lenkey^, and F. Kun^^ 

^ Department of Theoretical Physics, University of Debrecen, P.O. Box 5, H-fOlO Debrecen, 

Hungary 

^ Bay Zoltdn Nonprofit Ltd. for Applied Research, H-1116 Budapest, Fehervdri ut 130. Hungary 

We study the statistical features of magnetic noise accompanying the dynamic fracture of steel samples during 
mixed type fracture, where the overall ductile character of crack propagation is interrupted by a sudden brittle 
jump. The structure of the voltage time series is investigated by identifying discrete peaks which correspond to 
elementary steps of the jerky cracking process. We show that the height, duration, area, and energy of peaks 
have power law distributions with exponents falling close to the corresponding values of pure ductile fracture. 
For single peaks a power law correlation of the height and area with the width is evidenced, however, the mixed 
nature of fracture gives rise to a crossover between two regimes of different exponents. The average pulse 
shape of micro-cracking events has a parabolic form with a right handed asymmetry similarly to quasi-static 
fracture propagation. The asymmetry emerges due to stress localization at the advancing crack front. 


PACS numbers: 46.50.-Fa, 64.60.av, 05.90.-Fm 

Dynamic fracture of heterogeneous materials plays a 
significant role in a broad range of industrial applications 
from machine engineering to the safety of large scale con¬ 
structions. For the assesment of the dynamic fracture 
strength of materials a detailed understanding of the ini¬ 
tiation and propagation of cracks is essential^. Contrary 
to static fracture, the precise experimental characteriza¬ 
tion of dynamic fracture toughness is still a challenging 
problem. One of the main issues of experimental tech¬ 
niques is the determination of the onset time of crack 
propagation which is typically done by crack-tip moni¬ 
toring using high speed photography, single wire fracture 
gages, or strain gagesi*^. Due to the disorder present at 
the micro- or meso-scales of construction materials the 
propagation of cracks has been found to be a jerky pro¬ 
cess composed of discrete jumps which give rise to the 
emission of acoustic waves^. Dynamic fracture propa¬ 
gation in ferromagnetic materials like steel is also accom¬ 
panied by magnetic emission (ME) which yields an im¬ 
portant alternative to acoustic testing^i^. Such crakling 
noise measurements can help to overcome the limitations 
of traditional techniques. Recently, the integrated ME 
signal has been suggested to be capable to detect frac¬ 
ture initiation, however, the reliability of the method is 
not satisfactory^!^. More information is encoded in the 
structure and evolution of the ME time series of crack¬ 
ling events which can complement studies with integrated 
signals^. To work out well-established experimental tech¬ 
niques based on ME analysis with a high reliability a 
comprehensive understanding of the statistics and dy¬ 
namic origin of ME signals is needed^. 

In the present paper we investigate dynamic crack 
propagation by means of the magnetic emission technique 
focusing on the case of mixed fracture where the overall 
ductile character of crack advancement is interrupted by 
a sudden jump characteristic for brittle fracture. We ex¬ 
tend the standard ME analysis of the integrated statis¬ 
tics of pulse quantities by studying the shape of single 
pulses and the correlations of pulse characteristics. We 
demonstrate that the mixed-mode fracture gives rise to a 
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FIG. 1. Force F and voltage V signal of magnetic emission 
measured during the impact fracture of a V-notched JRQ 
RPV specimen. The overall ductile propagation of the crack 
is interrupted by a short brittle period indicated by the sud¬ 
den force drop accompanied by an intensive generation of high 
voltage pulses. 


crossover in the width-height or width-area correlation of 
crackling pulses, i.e. power law correlations are pointed 
out, however, different exponents are obtained for small- 
and large-sized pulses, typical for ductile and brittle crack 
propagation, respectively. The average temporal profile 
of single bursts proved to have a parabolic shape with 
a right handed asymmetry which is the consequence of 
stress localization at the propagating crack front. 

In order to investigate dynamic fracture of steel, exper¬ 
iments were carried out by means of the Charpy impact 
machine^^— . The machine itself is a large size pendu¬ 
lum which is raised to a certain hight and hits the spec¬ 
imen at the bottom of its swing. The specimens had 
the dimensions 1cm x 1cm x 5cm according to the stan¬ 
dards of Charpy tests^. Steel specimens were carefully 
manufactured in order to minimize the initial mechani¬ 
cal stress and damage inside the material. The velocity 
of the hammer of the impact machine had the same value 
vq = 5.5m/s in all the experiments. In order to control 
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FIG. 2. Identification of peaks and the definition of peak 
quantities such as height H, duration St, and area A. 


the position of the crack initiated by the energetic hit the 
specimens were either pre-cracked or V-notched in the 
middle. Samples were made of two different types of steel, 
i.e. reactor pressure vessel steel JRQ RPV was used for 
both pre-cracked and V-notched specimens, while carbon 
steel S355J2 was only used for V-notched samples. Since 
the material of the samples has ferromagnetic properties 
the propagation of the crack is accompanied by changes 
of the magnetic field. The ME activity was recorded 
in the form of a voltage signal by a small-sized coil at¬ 
tached to the impactor. Recently, it has been shown that 
this magnetic noise provides interesting insight into the 
dynamics of crack propagation being sensitive to fluctua¬ 
tions of the speed of the crack front and of the opening of 
the cracb^. Parallel to magnetic noise we also measured 
the force acting on the moving impactor. The sampling 
rate was At = 5 x 10“^s which provided sufficient reso¬ 
lution at our impact velocity vq. 

Figure [T] presents an example of the voltage signal V (t) 
and of the force F{t) measured simultaneously during an 
impact fracture experiment. Note that the type of frac¬ 
ture, i.e. either brittle or ductile can be inferred from the 
F{t) curve: brittle fracture implies the sudden initiation 
of a crack running through the sample which gives rise 
to a sudden drop-down of the force after a short upload¬ 
ing. The gradual evolution of F{t) observed in the figure 
is characteristic for ductile fracture where stable crack 
propagation emerges. In our experiments all specimens 
had mixed-fracture, i.e. the crack propagation started in 
a ductile manner, however, the sudden fall of the force 
indicates a brittle jump of the crack which was then again 
followed by stable propagation (see Fig. [T]). In spite of 
the complex cracking process, the force curve varies rel¬ 
atively smoothly, the wavy pattern is caused by elastic 
waves traveling through the sample. However, the volt¬ 
age signal shows strong fluctuations with a complicated 
structure along the entire fracture process. Note that the 
brittle regime of fracture typically gives rise to a more in¬ 
tense magnetic activity with higher and broader peaks^. 


In order to characterize the structure of magnetic noise 
accompanying dynamic fracture we decompose the volt¬ 
age time series into a sequence of discrete peaks following 
the method presented in Ref^^. Note that the method 
is analogous to the evaluation technique of Barkhausen 
noise measurements^^— . As a first step we identify the 
background noise level 14 by searching for the most prob¬ 
able voltage value. Then peaks are identified as segments 
of the V (t) time series between two intersections with the 
background 14. In a single sample a few hundred volt¬ 
age peaks are identified where peak i starts at time if 
and ends at tf. In order to characterize single peaks a 
set of quantities are determined by carefully evaluating 
the measured voltage signals, see Fig. [2] for illustration: 
The overall geometry of peaks is characterized by their 
height Fli and width, i.e. the peak duration SU. The 
height Fli of peak i is the largest voltage value = 

max [V{t){tf < t < tf)] in the peak duration from which 
the background is subtracted = 14 "*“^ — 14 . The peak 
duration Sti is obtained as the difference of the ending 
and starting times if as Sti =lt~ti- The area Ai of peak 
i is calculated as the integral of the voltage curve above 
the background level (see Fig. HD, while for the energy of 
the peak Ei the square of the voltage signal is integrated 
in a similar way 

A^ = [' {V{t)-Vb)dt, E, = [' {V{t)-Vbfdt. (1) 

Jti Jt‘ 

It has been shown in Refi^ that peaks of magnetic noise 
can be associated to elementary events of dynamic crack 
propagation. Due to the stochastic nature of crack prop¬ 
agation all quantities defined above have strong fluctua¬ 
tions, hence, for the characterization of the voltage time 
series we determined the probability distributions and 
the correlations of peak quantities. In order to improve 
the statistics of the data experiments were repeated 5, 7, 
and 12 times for pre-cracked, V-notched JRQ RPV, and 
for V-notched carbon steel, respectively. 

The primary quantities characterizing the shape of 
peaks are the width St and height H. Figure[3Ka) presents 
the probability density of peak heights p{F{) for all three 
sets of experiments. Power law decay of p(H) can be 
observed for large values of H 

p{H) cx H-P, (2) 

where the exponent j3 has nearly the same value P = 
1.55 ± 0.15 for all types of samples. The corresponding 
distribution of peak durations p{St) has also a power law 
asymptotics 

p{St) oc St~°‘ (3) 

persisting over three decades in St (see Fig. [3K6)). The 
value of the exponent can be determined with a rela¬ 
tively good precision a = 1.67 ± 0.06. Note that due to 
the low statistics of the data in the large St regime it is 
hard to infer the shape of the cutoffs limiting the power 
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FIG. 3. Probability distributions of the height H (a), duration 
St (fe), area A (c), and energy E (d) of voltage peaks. Power 
law functional form is evidenced in all cases. 


laws. The area A and energy E defined by Eqs. ([T|) pro¬ 
vide measures for the magnitude of crackling events rep¬ 
resented by the peaks. The probability distributions of 
peak area p{A) and energy p{E) are presented in Figs. 
|3^c, d), where a power law functional form is evidenced 
for both distributions over several orders of magnitude 

p{A) oc A~'^, p{E) oc E~'^. (4) 

The exponents were obtained numerically as r = 1.3 ± 
0.05 and e = 1.2±0.05. The power law character of single 
peak quantities implies that the jerky steps of crack prop¬ 
agation are not random but correlations occur between 
consecutive jumps of the crack front. Ductile cracks ad¬ 
vance in small steps as the hammer of the impact machine 
proceeds. In case of brittle fracture the crack decouples 
from the impactor and it leads to sudden failure with a 
small number of large jumps^. Consequently, brittle and 
ductile cracking are characterized by lower and higher ex¬ 
ponents of the distributions of peak quantities^. A very 
interesting outcome of our measurements is that mixed 
type fracture is also characterized by power law distribu¬ 
tions which are naturally dominated by the peaks gener¬ 
ated in the ductile regime. 

Of course, single peak quantities are not independent 
of each other, i.e. it can be expected that a broader peak 
has a larger height and in turn a larger magnitude repre¬ 
sented by the peak area and energy. Such correlations of 
the peak height E[ and duration St can be quantified by 
the conditional average, i.e. averaging El for fixed values 
of St. Figure m^a) shows that the average height {H) has 
a power law dependence on the width of peaks St 

( 5 ) 




FIG. 4. Gorrelation of peak quantities. Average height {H) 
(a) and average area (A) (b) of peaks as a function of the 
peak duration St. Power law dependence is obtained with a 
crossover between two regimes of different exponents. 


however, the value of the exponent ^ has a distinct change 
at a characteristic peak duration St^. The crossover be¬ 
tween the two power law regimes of different exponents 
occurs at Stc ^ 2 x 10“®s and the exponent ^ takes the 
values ^ = 1.35 and ^ = 0.35 below and above Stc, respec¬ 
tively. The statistics of the data in Fig.lH^a) is somewhat 
lower above the crossover point St > Stc which is consis¬ 
tent with the rapid decrease of p{St) in Fig. [S^fo). The 
result implies that for broad peaks St > Stc the height in¬ 
creases slower with the width than for the narrow ones. 
Note that the crossover is not visible on the probabil¬ 
ity distributions p{H) and p{St). A possible explanation 
is that both the height and duration values above the 
crossover fall in the cutoff regime of the distributions. 

Based on Fq. we can assume the approximate re¬ 
lation H oc St^ from which a scaling relation of the expo¬ 
nents a and (3 of p{St) and p[H) can be derived as 

a = l + ^(^-l). (6) 

Substituting the value of the exponents the scaling form 
is fulfilled with a good precision below the crossover du¬ 
ration but not above Stc since the exponents a and /3 are 
dominated by the ductile regime of crack propagation. 
Based on the height-width relation of pulses the area A 
can be simply approximated as a function of St 

A oc St^~^^. (7) 

To test the validity of the above form we determined 
the conditional average (A) (St) by directly averaging A 
in bins of St. Figure \Mb) demonstrates that power law 
behavior is obtained as expected with the same crossover 
behavior that has been found for the peak height (see 
Fig-HKa))- The values of the exponent 2.3 and 1.3 below 
and above the crossover point are consistent with Eq. ©■ 

The correlation of the height, area, and duration pro¬ 
vides a characterization of the overall shape of single 
cracking peaks. A more detailed description can be 
obtained by determining the average temporal profile 
{V{t, St)) of pulses of a fixed duration S t^^3'^ . Since sin¬ 
gle voltage pulses correspond to individual steps of the 


(H) oc St^, 
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FIG. 5. Average temporal profile of pulses rescaled with N{St) 
varying the duration St/At in a broad range from 10 to 270. 
The continuous line represents best fit with Eq. 


advancing crack, the average temporal profile provides in¬ 
formation about the time evolution of intermittent crack¬ 
ling bursts. To improve the statistics, pulses of duration 
St ± 0.l6t were averaged together using the data of all 
measurements. Figure [5] presents that rescaling {V{t, 6t)) 
with the norm N{St) = {V(t/St)) d{t/6t), pulse pro¬ 
files of different duration St can be collapsed on a master 
curve as a function of the normalized time t/St. It follows 
from Eq. © that the norm scales as N{St) ~ St^ so that 
the good quality data collapse implies the validity of the 
scaling form of pulse shapes 

{V{t,St)) =St^ fit/St). ( 8 ) 

The scaling function f{x) has a parabolic shape which 
shows that during crackling jumps the crack front starts 
slowly then it accelerates followed by a gradual slow 
downi^iii. Note that the profile f{x) is slightly skewed 
to the right indicating the time irreversible nature of the 
temporal burst dynamics. The same temporal profile 
has been found very recently for crackling bursts gener¬ 
ated at a slowly propagating crack front in heterogeneous 
material a^^d^ . It has been shown that the symmetry of 
the pulse shape carries information about the range of 
stress redistribution following local breaking events: long 
range interaction (mean field) results in symmetric pro¬ 
files, while short range load sharing gives rise to a right 
handed asymmetryi^d^. It follows that the right handed 
asymmetry of the pulse shapes in Fig.[S]is the fingerprint 
of the localized stress redistribution along the crack front 
during dynamic fracture propagation. The scaling func¬ 
tion f(t/St) was fitted with the functional form 

f{x)^[x{l-x)r[l-a{x-l/2)], (9) 

where the exponent 7 determines the shape of f{x) in the 
vicinity of a: = 0 and a: = I, while a controls the degree of 


asymmetryi^. Best fit was obtained with 7 = 0.55 which 
is significantly lower than the corresponding mean field 
value 7 = however, it falls closer to the result 

on quasi-static planar crack propagation 7 Ri 0.78^^. The 
asymmetry is quantified by the parameter value a = — 0.2 
which has a larger absolute value than the ones predicted 
by models of static crack propagatio n^'^d^ . 

We studied dynamic fracture of steel by analysing the 
structure of magnetic emission time series generated by 
the jerky propagation and opening of a crack. Our anal¬ 
ysis revealed that for all types of specimen the height, 
duration, area, and energy of crackling pulses are all 
power law distributed with exponents falling close to 
their counterparts of pure ductile fracture. The mixed 
(ductile-brittle) nature of the fracture process gives rise 
to a crossover in the correlation of peak quantities with 
a sharp transition between two regimes of different ex¬ 
ponents. The most remarkable result is that the average 
temporal profile of crackling peaks has a parabolic form 
with a right handed asymmetry similarly to pulse shapes 
characterizing quasi-static crack propagation in hetero¬ 
geneous materials. Profiles of different durations can be 
collapsed on a universal master curve on both sides of the 
crossover point. The asymmetry of pulse shapes can be 
attributed to stress localization at the propagating front. 
The results imply that the time evolution of ME peaks 
is the same at the onset of crack propagation as later 
in the regime of stable propagation since it is solely de¬ 
termined by the range of stress redistribution. To obtain 
information on the onset time of fracture the analysis has 
to go beyond single peak quantities and involve e.g. the 
waiting time between consecutive crackling events. 
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